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Origin of strong magnetic fields in Milky-Way like galactic haloes 
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ABSTRACT 

An analytical model predicting the growth rates, the absolute growth times and the saturation 
values of the magnetic field strength within galactic haloes is presented. The analytical results 
are compared to cosmological MHD simulations of Milky-Way like galactic halo formation 
performed with the N-body / Spmhd code GADGET. The halo has a mass of fts 3 ■ 10 12 
Mq and a virial radius of « 270 kpc. The simulations in a ACDM cosmology also include 
radiative cooling, star formation, supernova feedback and the description of non-ideal MHD. 
A primordial magnetic seed field ranging from 10 -10 to 10 -34 G in strength agglomerates 
together with the gas within filaments and protohaloes. There, it is amplified within a couple 
of hundred million years up to equipartition with the corresponding turbulent energy. The 
magnetic field strength increases by turbulent small-scale dynamo action. The turbulence is 
generated by the gravitational collapse and by supernova feedback. Subsequently, a series of 
halo mergers leads to shock waves and amplification processes magnetizing the surrounding 
gas within a few billion years. At first, the magnetic energy grows on small scales and then 
self-organizes to larger scales. Magnetic field strengths of « 10~ 6 G are reached in the center 
of the halo and drop to ps 10~ 9 G in the IGM. 

Analyzing the saturation levels and growth rates, the model is able to describe the process of 
magnetic amplification notably well and confirms the results of the simulations. 

Key words: cosmology: early universe - galaxies: formation - galaxies: halos - galaxies: 
magnetic fields - methods: analytical - methods: numerical 



1 INTRODUCTION 

The Lambda Cold Dark Matter model (A CDM) is the standar d 
tool describing the evolution of the universe (Ko matsu et alj20lH) . 
Quantum fluctuations in the primordial energy distribution develop 
into the condensing matter and trigger gravitational instabilities. 
Dark matter clumps in filaments and protohaloes, and subsequently 
baryonic matter falls into the potential wells of the dark matter, 
thereby forming the first stars and galaxies. In a hierarchical pro- 
cess of merger events , larger structures grow dWhite & Reesll 19781 : 
White &Frenflll99lh . Numerical simulations of structure forma- 
tion within a ACDM universe show good agreement between the 
calculated and the observed distribution of matter and structures 
dSpringel et alj2005bl . ] 2006 ) . However, cosmic magnetic fields are 
still widely neglected in these kind of simulations, although their 
presence can influence the dynamics of an astrophysical system 
significantly. 

Observations reveal strong magnetic fields of /nG stren gth in late- 
type g a laxies (for reviews on cosmic magnet i sm se e e.g. lBeck et al.l 
( 1996), Widrow (2002), Ku lsrud & Zweibell d2008h . and references 
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therein). The energy density of these magnetic fields is compara- 
ble to other dynamically important energy densi ties, i.e. the mag- 
netic fi eld seems to be in equipartition with them. lNeronov & Vovkl 
(2010) also find strong magnetic fields permeating the intergalac- 
tic medium (IGM). The IGM magnetic fields are highly turbu- 
lent dRvu et all 20081) and th eir strength is estimated to the order 
of nG ( Kronberg et al. 2008). Additionally, mag netic fields of fj.G 
strength can be found in high redshift galaxies dBernet et alj2008h . 
Also, there is evidence of highly magnetized damped Lyman alpha 
systems at redshift m 2, w hich act as building blocks for galactic 
systems dWolfe et all2005h . 

The origin of these magnetic fields is still unclear. Global pri- 
mordial magnetic fields can be seeded by battery processes in 
the early universe dBie rmann 1950; Mishustin & Ruzmaikinl l 1 9721 : 
IZel'dovich et all 1 19831 : iHuba & Fedderlll993l) . Alternatively, seed 
fields can be generated by ph ase transitions after the BigBang or 
various other mechanisms (seelwidrow 2002, for a review). 
In a subsequent process, these weak seed fields of sometimes 
1Q~ 20 G have to be amplified to the observed values. 
ILesch & Chibal |~1995) demonstrate the possibility of strong mag- 
netic fields at high redshifts through battery processes and pro- 
togalactic shear flow amplification. The presence of strong /j,G 
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galactic magnetic fields is commonly explained by galactic dy- 
namos converting angular momentum into magnetic energy in dif- 
ferenti ally rotating disks. Th e two main theories are the a-uj dy- 
namo i Ruzmaikin et alj|l979[) or the cosmi c ray driven dynamo 
jLesch & Hanasj 120031: lHanasz et al] 120091). For re vie ws of dy- 
namo theory see lBrandenburg & Subramaniar] d2005h or lShukurovl 
j2007t) . However, these dynamos operate on timescales (e-folding 
time, not absolute amplification time) of the order of 10 8 yrs and re- 
quire a differentially rotating galactic disk. Hence, irregular galax- 
ies at high redshift have to be magnetized by another process. An- 
other po ssible magnetization proc ess is the "Cosmic dynamo" as 
given bv lDubois & Tevssierl d2010t) . Within their approach, the uni- 
verse is magnetized by gravitational instabilies, galactic dynamos 
and wind-driven outflows of gas at times of violent star formation 
activity. 

Small-scale dynamos operate on timescales of the order of 
10 6 yrs through random and turbulent shear flow motions 
jBiermann & Schliited 1 195 lh . Magnetic energy increases ex- 
ponentially on small scales first by stretching, twisting and 
folding the magnetic field lines by random motions and 
then organizing them on the largest turbulent eddy scale 
jZel'dovich et alj|l983l; iKulsrud & Andersorj|l992l: I Kulsrud ct al. 
19971; iMalvshkin & Kulsrudl2002l ; ISchekochihin et al.l2002Ll2004l ; 



Schleicher etai] T2Q10). Galaxy mergers are a natural part of the 



bottom-up pi c ture of the grow th of struc t ures i n the universe. 
iKotarba et af] d20ld 1201 lh and iGeng et ail ( 120121) show that tur- 
bulence induced during galactic major and minor mergers is able to 
amplify magnetic fields in galaxies and in the IGM up to equipar- 
tition between the magnetic and turbulent energy density, as ex- 
pected from the small-scale dynamo theory. This theory is a good 
method to desc ribe the amplification processes and corresponding 
timescales (e.g. Arsh akian et alj2 009). However, galactic dynamos 
are still inevitable to explain the regularity of galactic magnetic 
fields and their spiral structure, which are revealed by observations. 
Analytical calculations and cosmological simulations of structure 
formation including the evolution of magnetic fields can give new 
insights in the physical processes of creating, amplifying and satu- 
rating magnetic fields in the universe on all kind of scales. 
In this work, an analytical model predicting the growth rates, the 
absolute growth times and the saturation values of the magnetic 
field strength within galactic haloes is presented. The analytical re- 
sults are compared to cosmological MHD simulations of Milky- 
Way like galactic halo formation including star formation and non- 
ideal MHD. It is shown that the analytical model and the cosmolog- 
ical simulations agree notably well for different initial, primordial 
magnetic seed fields spanning a range of 25 orders of magnitudes. 
The paper is organized as follows: The analytical calculations 
are shown in section 2. Section 3 briefly describes the numerical 
method. In section 4 the cosmological inital conditions and the 
magnetic seed field are motivated. A detailed analysis of the per- 
formed simulations and the magnetic field amplification is given in 
section 5. Section 6 compares the numerical results with the ana- 
lytical description. The main results are summarized in section 7. 



hydrodynamical Reynolds numbers, a stationary flow transits from 
the laminar regime into the turbulent regime and becomes unstable. 
Hence, an overview of the characteristics of a non-stationary per- 
turbated magnetic field in such an unstable flow in a cosmological 
context is given. 



2.1 Local perturbation ansatz 

The Reynolds number is a characteristic dimensionless quantity de- 
scribing the ratio of inertial forces and viscous forces of a flow: 



Re = 



inertial forces 
viscous forces 



VL 

1 

V 



(1) 



with V and L being typical velocity and length scales of the flow 
and v the kinematic viscosity. A stationary flow will become un- 
stable if the Reynolds number Re exceeds a critical value i?e cr it- 
In this case, any initial infinitesimal perturbation will be amplified 
through the flow. In the following calculations, the hydrodynami- 
cal flow is assumed to be unstable, since the inertial length scale 
is sufficiently larger than the viscous length scale. The magnetic 
Reynolds number Rm, which is defined as Rm — VL/rj (with the 
magnetic resistivity rf), is also sufficiently large to allow for pertur- 
bations to grow. The induction equation of ideal MHD reads: 



dB(x,t) 
dt 



[v(x,t) x B(x,t)] . 



(2) 



Within the small perturbation approximation, the velocity field 
v(x,t) and the magnetic field B(x,t), respectively, can be de- 
composed into the sum of a stationary component Vq(x, t) and 
Bq(x, t), and a perturbated component vi(x, t) and Bi(x, t), re- 
spectively. For the induction equation (O, this procedure results in 
the following decomposition: 



dB 

dt 



— V x [vq x B ] 



dt 



V x [v x Bi + v± x Bo] 



(3) 



(4) 



The right crossproduct in Eq. $4$ can be dropped, because only the 
growth of magnetic perturbations in a stationary flow is important, 
and the growth of velocity perturbations in a weak stationary mag- 
netic field can be neglected. The general solution of Eq. ((4) is a 
sum of special solutions, whereby Bi includes a time-dependent 
factor e~ lSlt . The complex frequency Q. is given by Q — oj + il\ 
with periodicity u and growth rate F. For growing perturbations 
r will be positive. The perturbated component Bi can be further 
decomposed into 



Bx{x,t) = B t (t)s(x) 



(5) 



2 ANALYTICAL DESCRIPTION 

This section gives an analytical approach describing the behaviour 
of the magnetic field strength during halo formation. In order to 
derive an analytical model, expressions for the cosmological de- 
cay, the exponential amplification process, the saturation and the 
relaxing decay of the magnetic field strength are needed. For large 



B t {t) = B t {t )e vt e-^\ (6) 

with a spatially dependent complex vector function s(x) and a 
complex scalar amplitude Bt(t). The time derivative of the square 
of the amplitude is: 

^l^=2TBUt). (7) 
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Eg. Q was also derive d for a kinematic dynamo by 



lKulsrud& AndersorJ dl992h and for a turbulent magnetized 
dynamo by Malvs hkin & Kulsrud ( t2002h using spectral anal- 
ysis of the growth of the magnetic energy. In the weak-field 
approximation, the flow is able to amplify the frozen-in magnetic 
perturbations, since the magnetic field is frozen into the velocity 
field. Thereby, any information about the original orientation of 
the magnetic seed field is lost. As soon as equipartition is reached, 
the back-reaction of the magnetic field on the velocity field has 
to be considered and the weak-field approximation breaks. Then, 
one can either solve the full Navier-Stokes equations including 
magnetic pressure and tension forces or model the b ack-reaction 
by truncating the growth rate F. Belva nin et al.l (l993) describe the 
latter "Equipartition dynamo" approach by expanding V in a power 
series and considering second-order terms to truncate the growth 
rate. These non-linear effects force the maximum amplitude to an 
equipartition value and the dynamo process saturates. Generally, 
a "Supraequipartition dynamo" should be expected, because 



compression of B 



(e.g. the effective adiabatic index 



is 7adi = 4/3 for isotropic compression) can still occur. This 
would result in the magnetic energy density first to shoot over the 
turbulent energy density (7 a di = 1). and then decay towards the 
saturation level. However, this effect is too small to be significant. 
Now, to account for the saturati on of the magnetic field at the 
equipartition level, T is amended dBelvanin et al.l 1993) as follows 



7 



1 



B- 



(8) 



and the equation for the growth of the magnetic field amplitude 
takes the form: 



dBfjt) 
dt 



27 



Bt(t) 
B 2 

^sat 



(9) 



This growth of the magnetic perturbations is an iterative process: 
Bi can only grow to the order of magnitude of Bo, then a new Bo 
and Bi have to be defined. Furthermore, since turbulence is the 
driver of the magnetic field amplification, equipartition between the 
magnetic and the turbulent energy density is a good approximation 
for the saturation magnetic field (no further growth will occur, once 
the magnetic field has saturated): 



B. 



sat _ 1 2 
>7T Z 



(10) 



with the density p and the turbulent velocity Vturb of the system. 
Together wi th the initial condition B 2 (t = 0) = Bq Eq. ((9} has 
the solution dLandau & Lifshitzlll959h : 



Bt(t) 



(id 



As a next step the growth rate 7 has to be determined. The ki- 
netic energy of the turbulence can be assumed to follow a one- 
dimensional Kolmogorov spectrum of the form 

/(A) ~ ^u rb fc- 5/3 , (12) 

where ^tnrh is the m ean square turbulent velocity. From Eq. {7} 
iKulsrud et alj d 19971) find for the growth rate 7 



2 7 « J ^^-dk « J y/kl{k~jdk, 



(13) 



with an eddy turnover rate kv^ « [kl(k)] 1 ^ 2 and with Vk being 
the typical eddy velocity on k scale. Integrating Eq. J 1 3t yields the 
growt h rate 7. Performing additional calculations including resis- 
tivity, [&lsmdet^] d 19971 ) find an advanced version of 7: 



7 = 2.050 



3/2 ,1/2 

"turb K turb 
1/2 
^turb 



(14) 



Eq. d 14b describes the turnover rate of the smallest turbulent eddy. 
Magnetic resistivity, which is essential for topological changes of 
the magnetic field lines through reconnection on small scales and 
transfer of magnetic energy into internal energy, and which is ef- 
fectively lowering the growth rate, is already included via r) tU rb- 
The dynamo process increases the magnetic energy on the dissi- 
pation scale first with a rate of 27 and saturates at a value com- 
parable to the turbulent energy on small scales first. The mag- 
neti c energy is then tr ansferred to larger scales with a rate of 
1 7 dKulsrud et alJI 199% by Lorentz forces unwrapping the folded 
field lines on small scales in an inverse cascade process. Finally, 
the magnetic field reaches saturation at a value comparable to the 
turbulent energy in the larg est eddy ( Kulsrud & Andersonl 1 19921 : 
iMalvshkin & Kulsrud I2OO2I) . The growth ra te given by Eq. dl4T > 
is t aken from the kinemat i c dyn amo theory dKulsrud e t al. 1997), 
but IMalvshkin & Kulsrucj {2002) find that the calculation leading 
to Eq. d 14b still holds for the magnetized turbulent dynamo in the 
weak-field approximation. 

These small-scale dynamos operate whenever turbulent and ran- 
dom motions and shear fl ows are stretching, twi s ting and folding 
the m agnetic fields lines dZel'dovich et al.|[l983l : ISchleicher et al.l 
l20ld) . Field lines come close on small scales first and hence, the 
magnetic energy increases at first on small scales. The amplification 
is dominated by the turbulent dynamo action and effects resulting 
from compression can be neglected within this model. 



2.2 Cosmological evolution and decay 

The next step towards an analytical model of the magnetic field 
amplification is to modify the growth equation ((9) such that the 
scalefactor a(t) describing the expansion of the universe during 
its evolution is accounted for. For an isotropic and stationary mag- 
netic field, the magnetic flux has to be conserved when space is ex- 
panding. From the first law of thermodynamics, the proportionality 
between the scalefactor a and the energy density of electromag- 
netic fields (or ultrarelativistic particles) e ~ B 2 , can be derived 
dLongaiilll998l) : 



e~a . (15) 

Any cosmological equation for the evolution of the magnetic field 
amplitude has thus to account for the proportionality Bt ~ a~ 2 . 
In cosmological simulations, the scalefacor a(t) instead of physi- 
cal time t is the natural integration variable. Given a flat universe 
without any radiation pres sure, the evolu tion equation for the scale- 
factor takes the form (e.g. Longair 1998h 



dtHo = da(— — + a 2 Q.A_) 
a 



(16) 



with Om and Q.a being the density contributions of matter and the 
cosmological constant, and Ho the present-day Hubble constant. 
The employed parameters of this ACDM cosmology (T able [Jl are 
close to the observed values (see Komatsu et alj|201lh . With the 
initial condition a(0) =0 Eq. d 1 6b has the solution: 
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Figure 1. Analytical growth curves for the magnetic field amplitude as a function of redshift for different numerical parameters. In each panel one parameter 
is varied, while the other parameters are held constant. Differences in the growth rates and saturation levels can be seen. 



ACDM Cosmology Parameters 



Matter density 
Dark energy density 
Total density 
Hubble constant 



Ho 



0.3 
0.7 
1.0 

70kms~ 



1 Mpc" 



Table 1. Parameters used for the ACDM cosmology. 



4(a) = 



asinh 




Combining equations ( lilt . J 1 4b . dl5b and J17| ( finally results in 



(17) 



B t (a) 



1 r 2-1 

— (47r/n? tulb ) + B 
or L 



2 g -2 7 t(a) 



(18) 



Eq. J 1 8b is an approximation for the growth of the magnetic field 
strength during the matter-dominated and later phases of the Uni- 
verse. 

In the beginning, the term a{t)~ 2 dominates and results in the 
cosmological dip at high redshifts. Then, the magnetic energy in- 
creases with e 27t(a ' during the gravitational collapse and due to 
star formation induced turbulence. This growth stops when equipar- 
tition is reached, since non-linear effects are truncating 7. Finally, 



when the system is relaxing and the turbulence is decaying, also the 
magnetic field will decay as a(t)~ 2 corresponding to a power-law 
decay with B t ~ t _4//3 . This decay is slightly stronger compa red 
to the B t ~ t~ 5/4 decay given by lLandau & Lifshit j ( fl95^) or 
lGeorgeifl992l) for the final stages of decaying kinematic turbulence 
(assuming that the magnetic energy decays with the same power- 
law as the turbulent energy, which is an approximation). The final 
stage of decay is reached by the time, when the Reynolds number 
becomes sufficiently small. This happens when the back-reaction 
of the magnetic field on the velocity field suppresses turbulent 
motions and also magnetic amplification and the system relaxes. 
ISubramanian et all j2006h found the power-law decay to set in al- 
ready for Reynolds numbers still as large as Re ~ 100 in galaxy 
clusters. 



2.3 Numerical parameters 

Eq. dl8b contains four free parameters: Vtmb, hurb, ?7turb and 
p. The numerical values for these parameters have to be taken 
from observations or extracted from numerical simulations. Table 
[2]shows the values used in this work. These values correspond to a 
timescale (e-folding time) of r = I/7 w 90 Myr and a saturation 
value of i3 sa t ~ 0.1/iG at redshift zero. Note, that these numerical 
parameters are constant mean estimates and do not reflect the time 
and spatial details of the simulations, but nevertheless are a good 
approximation. 

Estimating values for «turb an d ^turb is quite challenging, as the 
densities, velocities and length scales in galactic haloes range over 
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Analytical Model Parameters 

Turbulent length 'turb 25 kpc 

Turbulent velocity ^turb 75 km s — 1 

Turbulent dissipation rfturb 10 30 cm 2 s — 1 

Gas density n g!lB 10 — 5 cm - 3 

Table 2. Parameters used for the calculation of the analytical growth func- 
tion (physical units). 



Multi-phase Model Parameters 



Gas consumption timescale isf 

Number density threshold nth 

Mass fraction of massive stars /3 

Evaporation parameter A 

Effective supernova temperature T$n 

Temperature of cold clouds Tec 



2.1 Gyr 
0.13 cm" 3 
10% 
1000 
10 s K 
1000 K 



Table 3. Parameters for the star formation model ISpringel & Hernquisj 
l2003al) used in the simulations. 



many orders of magnitude. Therefore, such values can only be as- 
sociated to typical values within galactic haloes (see section[6]for a 
discussion which values reflect best the simulations). 
Fig. Q] shows growth curves as a function of redshift for a wide 
range of values of these parameters. In each panel, three parame- 
ter are held constant (see Table |2), while the fourth parameter is 
varied. As indicated by the presented calculations, «turbi ^turb and 
?7turb affect the growth rate 7, while wturb and p affect the satu- 
ration value .Bsat. Depending on the parameter configuration, the 
timescale ranges from order of 10 6 yrs to 10 s yrs and the satura- 
tion value for the magnetic field strength varies from 1CP 5 G to 
10" 8 G. 



3 NUMERICAL METHODS 

The simulations in this work are performed with the Em- 
body / SPMHD code GADGET i Springel etall l2001al ; ISpringell 
l2005al ; iDolag & Stasvszvd 120091) . Gadget uses a formula- 
tion of Sph, in which both energy and entropy are conserved 



dSpringel & Her nquis t 2002). For recent reviews on the SPH and 
SPMHD methods, see[Springel (2010h andlPricd (2012). Addition- 



ally, SUB FIND dSpringel et alJl2001bljDolag et al.ll2009h is applied 
to identify haloes and subhaloes and to calculate their respective 
center locations and virial radii. 

A detailed description of the SPMHD impleme ntation of MHD and 
its ext ension to non ideal MHD can be found in lDolag & Stasvszvnl 
d2009h and iBonafede et al] feoill) . Here, the standard (direct) 
SPMHD implementation is used, where the induction equation 



dB 
~dt 



V x (v x B) + r/AB 



(19) 



is evolving the magnetic field, and al so a spatially constant mag- 
netic resistivity r\ is applied. Following Bon afede et al.l J201 lh . non- 
ideal resistivity is assumed to be driven by the turbulence within 



10 d 



which is 



the gas and 77 tur b is of the order of 
consistent with models of the central region s of galaxy clusters 
( Schlickeiser et al. 19871; iRebusco et al]|2006h . Therefore, the con- 
stant turbulent resistivity describes the magnetic field decay on sub- 
resolution scales and is many order of magnitudes larger compared 
to numerical or ohmic resistivity, which are hence not of interest 
for this work. 

The magnetic field back-reacts on the velocity field via the Lorentz 
force . To account for the tensile instability (see lDolag & Stasvszvnl 
120091 : |Pricell2012l . for details) in SPMHD, the unphysical numeri- 
cal divergence force is substracted from th e equation o f motion fol- 
lowing an approach bv lB0rve et al. I feOOll) . Similar to lKotarba et al.l 
(2010), a limiter is applied to ensure that the correction force does 
not exceed the Lorentz force to avoid instabilities. 
To ensure a proper evolution of the magnetic field in numerical sim- 



ulations, it is of fundamental interest to maintain the V ■ B con- 
straint. In particular, an erroreous calculation can lead to unphysi- 
cal sources and sinks of magnetic energy. The MHD GADGET code 
keeps th ese numerical errors to a m inimum. For a detailed discus- 
sion, see lDolag & Stasvszvnl ( 2009) and section [531 
The implementation of MHD in GADGET was successfully em- 
ployed for the study of the m agnetic field evolution during star for- 
mation (Biirzle et al. 2011a. b), in isolated dKotarba et al.l2009l) and 

^ 11 — T] — k 

interacting galaxies (Kotarba et al. 2010, 2011; Geng et al. 2012), 
and in g alaxy clusters ijD onnert et al. 2 0091) . 

Also, the lSpringel & Hemquistl ( l2003al) star formation model is ap- 
plied. It describes radiative cooling, UV background heating and 
supernova feedback in a consistent two-phase sub-resolution model 
for the interstellar medium. Cold clouds with a fixed temperature of 
Tec are embedded into a hot ambient medium at pressure equilib- 
rium. These cold clouds are evaporating with an efficiency param- 
eter of A and form stars on a timescale of tsF, once they reach 
a density threshold of p t h- A fraction j3 of these stars is expected 
to die instantly as supernovae, heating the gas with a temperature 
of Tsn- Additionally, the hot phase is loosing energy via cooling, 
which is modeled assuming a primordial gas composition (Hydro- 
gen 76% and Helium 24%) with a temperature floor of 50 K (for 
details see lKatz et alj |T996). The cooling only depends on density 
and temperature, but not on metallicity. This star formation model 
leads to a self-regulated cycle of cooling, star formation and feed- 
back in the gas. 

Table[3]shows the numerical values of these parameters used in the 
simulations, which are performed without galactic winds. These 
numbers are choosen to reproduce the Kennicutt-Schmid t law be- 
tween surface density and surface star formation rate dSchmidtl 
ll959l;lKenmcuttlll998f) . 

However, for simulations of the turbulent small-scale dynamo the 
precise details of the star formation scheme are largely unimpor- 
tant. Cooling is required to obtain higher gas densities and smaller 
spatial scales in order to start efficient dynamo action. Furthermore, 
feedback-driven turbulence will contribute to the gravitationally- 
driven turbulence and raise the growth rates of the magnetic field 
strength. 



4 SETUP 

4.1 Dark matter initial conditions 

The presented simu lations start f rom c osmological initial condi- 
tions introduced by IStoehr et ai] d2002h . Starting point is a large 
ACDM dark matter only simulation box run with the GADGET 
code at different resolutions. The index of the power spectrum of 
the initial fluctuations is n = 1 and the fluctuation amplitude is 
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Simulation setup 

Scenario SF/Cool. N Gas N DM M Gas M DM B°£?f ati ° n B| 



ga0.bxO 


yes 


68323 


68323 


2.6 


gaO.bxlO 


yes 


68323 


68323 


2.6 


ga0_bxl4 


yes 


68323 


68323 


2.6 


ga0.bxl8 


yes 


68323 


68323 


2.6 


ga0.byl8 


yes 


68323 


68323 


2.6 


ga0_bx22 


yes 


68323 


68323 


2.6 


ga0.bx26 


yes 


68323 


68323 


2.6 


ga0.bx30 


yes 


68323 


68323 


2.6 


ga0.bx34 


yes 


68323 


68323 


2.6 


gal.bxO 


yes 


637966 


637966 


2.8 


gal.bxl8 


yes 


637966 


637966 


2.8 


gal_bxl8_nosf 


no 


637966 


637966 


2.8 


ga2_bx0 


yes 


5953033 


5953033 


3.0 


ga2.bx!8 


yes 


5953033 


5953033 


3.0 



1O 7 M 1.4-10 8 M Q - OG 

10 7 Mq 1.4-10 8 Mq x 10~ 10 G 

10 7 M o 1.4-10 8 M q x 10~ 14 G 

10 7 M Q 1.4-10 8 M q x 10~ 18 G 

10 7 M 1.4-10 8 M Q y 10~ 18 G 

1O 7 M 1.4-10 8 M q x 10~ 22 G 

10 7 M Q 1.4-10 8 M q x 10" 26 G 

10 7 M 1.4 ■ 10 8 M Q x 10~ 30 G 

10 7 M Q 1.4 ■ 10 8 M Q x 10~ 34 G 



10 6 M 1.5-10 7 M q - OG 
10 6 M 1.5 ■ 10 7 M Q x 10~ 18 G 

10 6 M Q 1.5 ■ 10 7 M Q x 10~ 18 G 



10 s M 1.6-10 6 M q - OG 
10 5 M Q 1.6-10 b M q x 10~ 18 G 



Table 4. Setup of the different simulations: The table lists wether star formation and cooling (SF/Cool.) is applied, the number of gas (Ncas) and dark matter 
particles (Nqm), the mass of the gas (M(j as ) and dark matter particles (Mtjm). as wel1 as the initial magnetic field orientation and strength for all simulated 
scenarios, respectively. 



erg = 0.9 (see lStoehr et alj |2002. for details). In a "typical" region 
of the universe a Milky-Way like dark matter halo is identified. The 
resulting simulations (with increasing resolution) are labelled GA0, 
GA1 and GA2, and contain 13 603, 123 775 and 1 055 083 dark 
matter particles, respectively, inside R200, which is the radius en- 
closing a mean density 200 times the critical density (virial radius). 
The forming dark matter halo is comparable to the halo of the 
Milky-Way in mass (« 3 x 10 12 M Q ) and size (» 270kpc). The 
halo was selected to have no major merger after a redshift of w 1 
and also a subhalo population comparable to the satellite population 
of the Milky- Way was fo und. More details about the prope r ties of 
this halo can be f ound in lStoehr etail d2002l) . Istoehr et alj ( 120031) 
and lStoehrl d2006h . Hence, GA0, GA1 and GA2 provide ideal ini- 
tial conditions to investigate the evolution of magnetic fields in a 
galactic halo, similar to the Milky- Way. 



4.2 Gas and magnetic field 

To add a baryonic component, the high resolution dark matter par- 
ticles are split into an equal amount of gas and DM particles. The 
mass of the initial DM particle is splitted according to the cosmic 
baryon fraction, conserving the center of mass and the momentum 
of the parent DM particle. The new particles are displaced by half 
the mean interparticle distance. 

The origin of cosmic magnetic fields is still unclear. During the 
early evolution of the universe, magnetic seed fields must have 
been generated by non-ideal mechanisms, which are independent 
of the magnetic field itself. It this work, primordial magnetic fields 
permeating the e ntire universe are assumed to be generated by 
battery processes teiermanrJll950l:lMishustin & R uzmaikini [l972l ; 

dovich et al. 1983; Huba & Fedder 1993) onset of structure for- 
mation. When deriving the induction Eq. l |19t of ideal MHD, the 
fluid is treated as a one particle-type plasma. Hence, in Ohm's law 
aE = j (with E the electric field and a the conductivity of the 
plasma), the current density j is only described by the motion of 
the protons. Strictly, the plasma also contains electrons and the cur- 
rent density j is a combination of the proton current density j 



and an electron current density j e . Then, electrons and protons are 
moving at different speeds in the plasma, reacting differently to 
perturbations. This results in currents and a non-ideal term in the 
induction equation of the form: 



dB(x,t) 
dt 



Vn c x Vp c 



(20) 



with the speed of light c, the elementary charge e, the electron den- 
sity n e , and the electron pressure p c . In shocks or non-isotrop re- 
gions, the gradients of n c and p c may be non-parallel and a mag- 



netic seed field of the strength 



10' 



G may be generated 



dBiermannlll950l; iMishusti n & Ruz maikinlll972l : Eel 'dovich et al.l 
1983; Huba & Fedder 1993). To seed the magnetic field within the 
simulations, a magnetic field vector is given to every gas particle 
pointing into the same direction carrying the same amplitude (most 
used choice: 10~ 18 G in x direction). This uniform setup gives an 
initially divergence-free magnetic field inside the simulation vol- 
ume. 

Actually, the magnetic energy should be distributed to the different 
scales of the simulation, resulting in a magnetic spectrum. How- 
ever, this spectral property of the magnetic field can be neglected 
for magnetic energy densites £ mag = B 2 /8n sufficiently smaller 
than the kinetic energy density £kin = pv 2 (weak-field approxi- 
mation). In this limit, the magnetic field B will be frozen into the 
velocity field v and follow its evolution. Hence, only the amplitude 
of the magnetic seed field is relevant, but not its direction. 
Table [4] shows the performed simulations, which can be grouped 
into three sets: First, low-resolution simulations GA0 are used for 
a numerical study of different seed field strengths ranging from 
10~ 10 G to 10" 34 G and being uniform in x direction. Addition- 
ally, a run of GA0 with a seed field in y direction is shown to con- 
firm the neglectability of the magnetic seed field direction. Second, 
for the seed field with the strength of 10 -18 G, higher resolution 
runs GA1 and GA2 are added to analyze the influence of the nu- 
merical resolution on the evolution of the magnetic field. Further- 
more, a special run of GA1 is performed to study the influence of 
the modeled star formation on the evolution of the magnetic field. 
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Finally, to study the influence of magnetic fields on the existing 
simulations, additional runs with only hydrodynamics are obtained. 



5 SIMULATIONS 

This section presents the results obtained from the numerical sim- 
ulations introduced above. Contour images of different quantities 
are created by projecting the SPMHD data in a comoving (1 Mpc) 3 
cube on a 512 2 grid using the code P-SMAC2 (Donnert et al., in 
preparation). 



5.1 Morphological and magnetic evolution 

An overview of the different stages of structure formation and 
their implications on the magnetic field is shown in Fig. [2] First, 
dark matter protohaloes and filaments are formed at redshift z ~ 
30 — 10, in which potential wells baryonic matter falls in. Within 
these structures the frozen-in magnetic field gets compressed. For 
isotropic compression this leads to B ~ p 2 ^ A . Note, that perpendic- 
ular compression of the magnetic field lines would lead to B ~ p. 
Due to the cosmological expansion, the magnetic field strength out- 
side these protohaloes decreases with a~ 2 . Turbulence is mainly 
created by the gravitational collapse. Secondly, as the gas density 
increases in the first structures, the threshold density is reached and 
star formation sets in at redshift « 10, further enhancing the ex- 
isting turbulence and consuming the available gas. In these dense 
regions, small-scale dynamo action starts, increasing the magnetic 
field strength exponentially, i.e. B ~ e 7 ' with the growth rate 7 
(see also section O. Merger events will create Shockwaves prop- 
agating into the IGM, possibly amplifying the magnetic field by 
compression and small-scale dynamo action. As equipartition is 
reached, the system relaxes and turbulent motions will decay, with 
additionally the magnetic field decaying. 

Fig. [H shows the projected number density of the gas at three dif- 
ferent redshifts in the simulation ga2_bxl8, together with the cor- 
responding virial radii of the forming galactic halo. Fig. [4] shows 
the corresponding projected total magnetic field strength, as well 
as arrows indicating the direction of the magnetic field. The differ- 
ent phases during the formation of the halo and the magnetic field 
amplification can clearly be seen. The magnetic field agglomerates 
together with the gas in filaments and protohaloes, where small- 
scale dynamo action is taking place. Furthermore, as merger events 
take place, Shockwaves are propagating into the IGM creating tur- 
bulence. The IGM magnetic field is amplified in stages with several 
Shockwaves propagating into it. Within each Shockwave, the mag- 
netic field is possibly amplified by compression within the shock- 
front and by small- scale dynamo action behind the shockfront (see 
iKotarba et ai]|2oT l], for an analysis of Shockwaves and their effect 
on the magnetic field during merger events). At redshift w 1 the last 
major merger event takes place and the magnetic field saturates, i.e. 
it evolves into energy density equipartition. The magnetic field sat- 
urates inside the halo at a mean value of w 10 -7 G and within the 
IGM at « 10~ 9 G, both at redshift 0. 

Fig. [3 shows the RMS magnetic field strength inside the galactic 
halo as a function of redshift for different seed field strengths. The 
amplification timescale (i.e. the gradient of the B(a) function dur- 
ing the exponential amplification phase) is the same for all seed 
fields and only the total time until saturation varies. After satura- 
tion, the magnetic field strength decreases with a power-law slope 
of « — 1, as irregularities in the magnetic field are dissipated. 
Fig. [6] shows radial profiles of the volume- weighted magnetic field 
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Figure 5. Growth curves of the volume-weighted RMS magnetic field 
strength inside the halo in the simulations GA0 for different seed fields. 
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Figure 6. Radial profiles of the mean magnetic field strength inside the halo 
in the simulation ga2_bxl8 for redshifts 1 and 0, respectively. For 2 = 1, 
which is just after the last major merger event, the field strength decreases 
with a slope of - 1 . 1 . For the relaxated system at 2 = the slope is -0.9. 



strength inside the galactic halo for two different redshifts. At red- 
shift 1 (after the last major merger), magnetic field strengths of 
several 10" 5 G are reached in the center of the halo and drop to 
~ 10~ 7 G at the virial radius with a slope of « — 1.1. At red- 
shift (virialized system with decaying turbulence), magnetic field 
strengths of several 10~ 6 G are reached in the center of the halo and 
drop to « 10~ 9 at the virial radius with a slope of w — 0.9. Since 
the gas density scales linearly with the distance from the galactic 
halo center, this indicates a relation of the form B ~ p. 
Summing up, within the presented simulations it is possible to am- 
plify a weak primordial magnetic field up to the observed equipar- 
tition values. 
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Figure 2. Overview of the different physical processes during structure formation. This table shows the stages of the evolution (column 1 ), during which an 
astrophysical process (column 2) is triggering a MHD mechanism (column 3) operating on the magnetic field, resulting in the equations and magnetic field 
strength values given in column 4. 
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Figure 3. Projected number density ra ga s in comoving units at different redshifts in the simulation ga2_bxl8. The shown regions are cubes with 1 Mpc 
(comoving) edge length centered on the halo center of mass. The white circles indicate the virial radius of the halo. The formation of filaments and protohaloes 
with subsequent merger events can be seen. 



5.2 Pressures and star formation 

Fig. [7] summarizes the energy flow and its effect on the star for- 
mation within the simulations. Via internal energy, star formation 
provides a sink (cooling) and source (SN injection) of internal en- 
ergy of the gas. Internal and kinetic energy are mutually exchang- 
ing via pressure forces and viscosity. Additionally, the gravitational 
collapse transforms potential energy into kinetic energy, which is 
converted partly back into potential energy through the fluid mo- 
tions in the potential wells. Furthermore, kinetic motions create 
magnetic energy (induction equation). The Lorentz-force describes 
the back-reaction of the magnetic field on the velocity field. Non- 
ideal resistivity redistributes the magnetic energy and also converts 
it into internal energy. Internal, kinetic and magnetic energy den- 
sities contribute to a total pressure. The balance between the grav- 
itational collapse and the total pressure support regulates the star 
formation rate. 



Fig. [8] shows the magnetic energy density 



= B 2 /8iv, the 



kinetic energy density e^m = P« 2 /2, the thermal energy density 
= (7 — l)pu and the turbulent energy density e tur b = 
P w turb/2 in the simulation ga2_bxl8. The adiabat ic index 7 is 
5/3 an d u denotes the internal energy. Similar as in lKotarba et al.l 
(2010), uturb is taken as an estimate of the turbulent velocity within 
the volume defined by a SPMHD particle. They find it to be a good 
SPMHD approximation of the turbulent velocity, although it over- 
estimates the turbulence on small scales and ignores the turbulence 
on scales larger than the smoothing scale. 

As shown in Fig. [8] the magnetic energy density increases from 
the seed value up to equipartition with the turbulent and thermal 
energy densities until a redshift of ~ 3. The cosmological dip can 
be clearly seen at the start of the simulations. The magnetic en- 
ergy density overshoots the turbulent energy density slightly in the 
beginning, which results from possible further compression after 
equipartition is reached. Afterwards, the virialized system relaxes, 
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Figure 4. Projected total magnetic field strength and magnetic field vectors in physical units at different redshifts in the simulation ga2_bxl8. The shown 
regions are cubes with 1 Mpc (comoving) edge length centered on the halo center of mass. Clumping of the magnetic field together with the gas in filaments 
and amplification within protohaloes can be seen. Furthermore, Shockwaves are driven into the IGM increasing the magnetic field strength, until it saturates 
on all scales. 



and the turbulent and magnetic energies decline. The thermal en- 
ergy density still rises, as the magnetic and turbulent energy are 
converted into thermal energy by resistivity and viscosity. With 
some delay, equipartition is also reached in the IGM by a stepwise 
amplification process: First, the magnetic field is amplified inside 
the most dense structures, and subsequently the IGM magnetic field 
undergoes merger-driven shock amplification. 

Fig. [9] shows the total star formation rate as a function of the cos- 
mological scalefactor for all GA0 simulations (see Tabled with 
different seed field strengths. Before the formation of filaments and 
protohaloes, no star formation takes place. At the time the first 
structures reach the necessary critical density (z w 10), star for- 
mation begins. As more gas is accreted onto the main halo through 
gravitational infall or due to merger events, the star formation rate 
rises. At a redshift of « 3 it peaks and then starts to decline. For 
the simulations without any magnetic fields, star formation is still 
constantly ongoing at a low rate by the end of the simulations. For 
the simulations including magnetic fields, the star formation rate 
decreases earlier than in the non-magnetized comparison runs, as 
soon as equipartition is reached (« fs 3), and is thus comperatively 
lower than in the comparison run or even stops completely by the 
end of the simulations. The magnetic configuration at equipartition 
provides additional support against further gas accretion, prevent- 
ing the gas inside the halo from reaching the threshold density nec- 
essary to form stars. 

Note, that the details of star formation depen d on the resolution of 
the simulation. ISpringel & Hemquistl J2003bh show that when in- 
creasing the resolution of the simulation, the time when the first 
stars form is pushed towards higher redshifts. Nevertheless, at low 
redshifts, the star formation rate converges to a value indepen- 
dent of the resolution. Furthermore, the total formed stellar mass 
does not change with resolution. These details in the star forma- 
tion model also influence the turbulent dynamo action. The start- 
ing point of the dynamo action depends strongly on the time when 
the first gas has collapsed, cooled and reached small enough spa- 
tial scales. As soon as star formation sets in, additional supernova 
energy is injected into the system, leading to more turbulence and 
higher magnetic field growth rates. The saturation value of the mag- 
netic field strength depends on the turbulent energy density, which 
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Figure 8. Volume-weighted energy densities as a function of redshift in the 
simulation ga2_bxl8 inside the halo and within the IGM. The magnetic en- 
ergy density (black line) gets amplified during the phase of halo formation 
until it reaches equipartition with the other energy densities, particularly the 
turbulent energy density (blue line). 



in turn depends on the the total injected feedback, and hence the 
total formed stellar mass. 

The simulation gal_bxl8_nosf (not shown) is performed with the 
same setup and methods, but with disabled star formation module. 
Within this simulation, the magnetic energy density does not get 
amplified up to equipartition with the turbulent energy density, but 
only rises a few orders of magnitude. This is clearly, because radia- 
tive cooling lowers the internal energy of the gas, thus allowing the 
gas to clump more heavily and reach higher densities. This results 
in smaller SPMHD smoothing lengths, and hence the small-scale 
dynamo will also operate on smaller scales, thus leading to higher 
growth rates. Summing up, radiative cooling and supernova feed- 
back are important in MHD simulations of galactic halo formation. 
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Figure 7. Diagram visualizing the flow of energy in the simulated MHD system. The star formation model provides a source (SN injection) and a sink (cooling) 
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Figure 9. Total star formation rate as a function of redshift in the simula- 
tions GA0 with different magnetic seed field. For simulations with magnetic 
fields, the star formation rate decreases when equipartition is reached and 
the additional magnetic pressure prevents the gas from reaching the density 
threshold required for star formation. 



5.3 Numerical reliability 

The numerical magnetic divergence (h | V • B\ /\B\) is a co mmon 
meas ure regarding the reliability of SPMHD simulations (e.g. lPricel 
l2012t) . It is calculated for every particle i inside its kernel, which 
is a sphere with radius equal to the smoothing length hi. Even for 
simulations employing the Euler potentials, which are free of phys- 
ical divergence by construction, this m easure can reach values of 
the order of unity dKotarba et ai1l2009h . The numerical divergence 



Figure 10. Mean numerical magnetic divergence measure 
(h | V • B\ /\B\) in ga2.bxl8 inside the halo and within the IGM. 



can be regarded as a measure for quality of the numerical calcula- 
tions and the irregularity of the magnetic field insi de each kernel 
and is not related to pos sible physical divergence iKotarba et al.l 
l20ld :l iBiirzle et alj|201 lal) . Here, an estimator for the numerical di- 
vergence of the form 

NumDivBi = V J 1 ' + % . ^ (B, - Bj) ■ VW{n h hi) 

t-r 1 \Bi\ + \Bj\ pi 

j 

(21) 

is used, where W is the SPMHD kernel function between the parti- 
cle i and its neighbours j. 
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Figure 11. Analytical growth functions (solid lines) as given by Eq. 1181 
using the parameters listed in Table [2] for different magnetic seed field 
strengths. After an initial cosmological dip, the magnetic field is grow- 
ing exponentially until it reaches equipartition with the turbulent energy 
density. Additionally, the simulated growth craves from Fig. [3] are shown 
(dashed lines). 



Figure 12. Volume-weighted RMS magnetic field strength inside the halo 
as a function of redshift for three simulations with the same magnetic 
seed field of 10~ 18 G but different resolution, together with the analyti- 
cal growth function (see Eq. U8t and Table[5J. The analytical function fits 
the simulated evolution of the halo magnetic field very well. 



Fig.[l0]shows the mean numerical divergence within the halo (solid 
line) and within the IGM (dashed line) as a function of scalefactor 
for the simulation ga2_bxl8. Throughout the entire simulation, the 
numerical divergence remains below unity. The numerical diver- 
gence is zero in the beginning of the simulations, as expected for 
a uniform magnetic seed field. During the phases of merger events 
and magnetic field amplification, the error estimator rises. The tur- 
bulent dynamo tangles the magnetic field lines and creates irregu- 
larities also below smoothing scales, resulting in a non-vanishing 
numerical divergence. During the phase of relaxation of the halo, 
NumDivB decreases, as the field lines are unfolded and ordered 
and the magnetic energy is dissipated or transferred to larger scales. 
Within the IGM, the numerical divergence remains constant. This 
is because within the IGM the NumDivB-decreasing process of 
magnetic field reordering is balanced by NumDivB-increasing pro- 
cesses. These are the accretion of gas onto the halo and the expan- 
sion of space, both resulting in an increasing smoothing length and 
thus an increasing NumDivB. 



6 AGREEMENT OF MODEL AND SIMULATIONS 

Since the simulations start at a finite redshift and not at 2 = 00, 
a simple shift of Bq — a"J.* rt j9s tart for the initial magnetic field 
strength in Eq. J 1 8b is used. The characteristic turbulent quantities 
in Eq. d 14b are assumed to be constant in space and time. This ap- 
proximation is good for the phase of strong star formation between 
redshift 10 and redshift 1, where also the majority of the magnetic 
amplification is taking place, while later the growth rate is trun- 
cated, and the form of the turbulent quantities is negligible. 
In Fig. QT] simulated growth curves of the magnetic field strength 
are shown (dashed lines). These curves match notably well with 
the calculated curves (solid lines). Table [2] shows the numerical 
values used for the calculations, which are resulting in a timescale 
(e-folding time) of w 90 Myr for the growth of the large scale mag- 
netic field. Extracting such values directly from the simulations is 
quite challenging, as the density, velocity and length scales within 



the simulated galactic haloes range over many orders of magnitude. 
However, the length scale l tU ib on which the magnetic energy den- 
sity increases first within the main halo is determined by the size 
of star forming regions, which (within the main halo) can be as- 
sociated to substructures. Such substructures (clum ps of gas, stars 
and dark matter) can be identified using SUBFIND dSpringel et al.l 
l2001bl : lDolag et al.l l2009). Only the largest substructures still con- 
tain gas and form stars within the main halo, having masses of 
« (10 s - 1O 1O )M and thereby diameters of a few 10 kpc. They 
are orbiting with typical velocities of a few 100 km s _1 , induc- 
ing gas RMS velocities within the main halo between 50 and 100 
km s _1 during the time of rapid growth of the halo, where also 
the magnetic amplification is taking place (between z ~ 10 and 
z ~ 1). Such values strongly motivate the choice of Zturb ~ 25 kpc 
and ftm-b ~ 75 km s _1 leading to the good agreement between the 
simulations and the analytical model. 

Fig. [12] shows a comparison of the analytical growth function (red 
line) as calculated according to Eq.[T8]using the values given in Ta- 
ble [2] together with the halo RMS magnetic field strength in the 
simulations ga0_bxl8, gal_bxl8 and ga2_bxl8, respectively. All 
curves start with the same magnetic seed field of B s t ar t = 10~ 18 
G. The resolution of the simulations increases by each a factor of 
roughly 10 from gaO to gal and from gal to ga2, respectively. 
Compared to the simulation with the standard resolution (brown 
line), the analytical growth function fits the simulated magnetic 
field evolution very well. Also, for higher resolutions (black and 
green lines), the fit is convincing. Note, that for higher resolutions, 
star formation sets in at higher redshifts leading to an earlier rise 
of collapsed, cooled gas and feedback and hence an earlier starting 
point for the turbulent dynamo. Nevertheless, the saturation value 
of the magnetic field strength is indistinguishable. 
For the calculation of the hydrodynamic and magnetic Reynolds 
numbers, the typical velocity V and length scale L are determined 
by the physical properties of the system, i.e. the soundspeed d? = 
7(7 — l)u and the Alfvenspeed v 2 — B 2 /flop, and the size of 
the halo, respectively, which do not depend on resolution. Within 
all the simulations, a constant turbulent resistivity r/ tur b is used, and 
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he nce, Rm stay s constant. Additionally, artificial viscosity as given 
bv lPricej d2012h is applied, where v depends on particle properties 
and spacings. A higher resolution thus leads to smaller v, and thus 
(given the constant V and L) to higher Re numbers and hence more 
turbulence. A higher (or, better resolved) turbulence in turn results 
in a higher growth rate of the magnetic field. 
The cosmological turbulent dynamo as described by Eq. l l 1 .St repro- 
duces the main features of the simulated non-ideal magnetic field 
amplification very well. 



7 SUMMARY 

In this paper the evolution of magnetic fields during galactic halo 
formation is discussed. The main focus is placed on the investi- 
gation of the processes responsible for the amplification of mag- 
netic fields from seed field levels to observed values. An analytical 
model for the evolution of the magnetic field driven by a turbulent 
dynamo is presented and the predictions of this model are com- 
pared with numerical, cosmological simulations of Milky-Way like 
galactic halo formation including the evolution of magnetic fields, 
radiative cooling and star formation. The most important results are 
summarized as follows: 

• A primordial magnetic seed field of low strength can be am- 
plified up to equipartition with other energy densities during the 
formation and virialization of a galactic halo in a ACDM universe. 
The final magnetic field strength decreases with a slope of « -1.0 
from w 10~ 6 G in the center to w 10~ 9 G behind the virial ra- 
dius (IGM) of the halo and also reaches as 10 G in interacting 
systems . These values are in notably good agr eement with obser- 
vations jBeck et al.lll996l:lKronberg et alj2008h . 

• The magnetic field amplification in filaments and protohaloes 
is dominated by turbulent dynamo action. Radiative cooling of 
the primordial gas is needed in order to reach spatial scales small 
enough for the turbulent dynamo to operate efficiently. Equiparti- 
tion is reached on small scales first and later on larger scales, con- 
sistent with theoretical expectations dBrandenburg & Subramanianl 
l2005l ; IShukurovll2007l : lArshakian et al.l |2009|> . The turbulence is 
driven by the gravitational collaps, by supernova activity and by 
mergers of protohaloes into the main galactic halo. After equipar- 
tition is reached, the magnetic energy decays with a power-law de- 
pendance of B t ~ t _4//3 . The IGM magnetic field is amplified by 
outflows of magnetized gas from the center of the haloes and by 
merger-driven shock amplification outside the main halo. 

• The amplification timescale (e-folding time) of the order of 
10 7 yr is small enough to describe the generation of strong mag- 
netic fields in irre gular galaxies at high redshifts as observed (e.g. 
iBernet et al.ll2008t> . 

• The structure of the resulting magnetic field is random and 
turbulent. Additional dynamo processe s, e.g. the a-uj dynamo 
dRuzmaikin et al.lll979l;lshukurovll2007h or the cosm ic-ray-driven 
dynamo dLesch & Hanaszll2003l ; lHanasz et al ] |200H are needed to 
produce regularity in the magnetic field topology. 

• Last, but not least, a basic analytical model is able to re- 
construct the numerical results very accurately. Weak magnetic 
perturbations grow in a non- stationary turbulent hydrodynamical 
flow. This amplification is slightly modified by the ACDM cosmol- 
ogy, particularly in the early universe. Non-ideal truncation of the 
growth rate finally yields equipartition. These processes result in an 
analytical growth function which fits the simulations astonishingly 
well. 



In the current picture of galaxy formation - together with cooling 
and star formation - magnetic fields are efficiently amplified from 
seed field levels up to the observed values. However, their detailed 
influence on the dynamics of the gas and the underlying seeding 
mechanisms still remain unclear and need to be investigated further. 
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